An Analytical and numerical investigation of escape rate for a 

noise driven bath 

Jyotipratim Ray Chaudhuri*, Suman Kumar Banik, Bidhan Chandra Bag and 

Deb Shankar Ray ^ 

Indian Association for the Cultivation of Science, Jadavpur, Calcutta 700 032, India. 

(February 1, 2008) 

Abstract 

We consider a system-reservoir model where the reservoir is modulated by 
an external noise. Both the internal noise of the reservoir and the external 
noise are stationary, Gaussian and are characterized by arbitrary decaying 
correlation functions. Based on a relation between the dissipation of the 
system and the response function of the reservoir driven by external noise we 
numerically examine the model using a full bistable potential to show that 
one can recover the turn-over features of the usual Kramers' dynamics when 
the external noise modulates the reservoir rather than the system directly. 
We derive the generalized Kramers' rate for this nonequilibrium open system. 
The theoretical results are verified by numerical simulation. 
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I. INTRODUCTION 



More than half a century ago Kramers proposed a diffusion model for chemical reactions 
in terms of the theory of Brownian motion in phase space Since then the model and 
several of its variants have been ubiquitous in many areas of physics, chemistry and biology 
for understanding the nature of activated processes in classical quantum and semiclas- 

sical PHTTH systems, in general. These have become the subject of several reviews [P"^-P^ 
and monograph ||T3| in the recent past. 

In the majority of these treatments one is essentially concerned with an equilibrium 
thermal bath at a finite temperature which stimulates the reaction co-ordinate to cross the 
activation energy barrier. The inherent noise of the medium is of internal origin. This im- 
plies that the dissipative force which the system experiences in course of its motion in the 
medium and the stochastic force acting on the system as a result of random impact from 
the constituents of the medium arise from a common mechanism. From a microscopic point 
of view the system-reservoir Hamiltonian description [p!^-[T9| developed over the decades 
suggests that the coupling of the system and the reservoir co-ordinates determines both 
the noise and the dissipative terms in the Langevin equation describing the motion of the 
system. It is therefore not difficult to anticipate that these two entities get related through 
a fluctuation-dissipation relation ( these systems are sometimes classified as thermody- 
namically closed system in contrast to the systems driven by external noise in nonequilibrium 



statistical mechanics ||2T| ). However, when the reservoir is modulated by an external noise 
it is likely that this relation gets affected in a significant way. Since the modulation of the 
reservoir crucially depend on its response function, one can further envisage a connection 
between the dissipation of the system and the response function of the reservoir due to the 
external noise from a microscopic point of view. 

In the present paper we explore this connection in the context of activated rate processes 
when the reservoir is modulated by an external noise. Specifically our object here is twofold: 
(i) to explore the role of reservoir response as a function of external noise on the system 
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dynamics in contrast to direct driving of the system by external noise, (ii) to calculate the 
generalized Kramers' rate for the steady state of this nonequilibrium open system taking 
full care of thermodynamic consistency condition. Both the internal and the external noises 
are Gaussian, stationary and are characterized by arbitrary decaying correlation functions. 
While the internal noise of the reservoir is thermal, the external noise may be of thermal 
or non-thermal type. We consider the stochastic motion to be spatial diffusion limited and 
calculate the rate of escape in the intermediate to strong damping regime. It is worth to 
mention that the externally generated nonequilibrium fluctuations can bias the Brownian 
motion of a particle in an anisotropic medium and may be used for designing molecular 
motors and pumps [^. We further mention that nonequilibrium, non-thermal systems 
have also been investigated phenomenologically by a number of workers in several other 
contexts, e.g., for examining the role of color noise in stationary probabilities properties 



of nonlinear systems nature of crossover p5[, rate of diffusion limited coagulation 



processes [Q, effect of monochromatic noise etc. While these treatments concern 
direct driving of the system by an external noise the present consideration is based on 
modulation of the bath. A number of different situations depicting the modulation of the 
bath by an external noise may be physically relevant. As, for example, we consider a simple 
unimolecular conversion (say, an isomerization reaction) from A — > B. The reaction is 
carried out in a photochemically active solvent under the influence of external fluctuating 
light intensity. Since the fluctuations in the light intensity result in the fluctuations in the 
polarization of the solvent molecules, the effective reaction field around the reactant system 
gets modified. Provided the required stationarity of this nonequilibrium open system is 
maintained (which is not difficult in view of the experiments performed in the studies of 
external noise- induced transitions in photochemical systems ||2^) the dynamics of barrier 
crossing becomes amenable to the present theoretical analysis that follows. 

The remaining part of this paper is organized as follows: In Sec. II we discuss a system- 
reservoir model where the later is modulated by an external noise and establish an interesting 
connection between the dissipation of the system and the response function of the reservoir 



due to external noise. The stochastic motion in a hnearized potential field is described in 
terms of a Fokker-Planck equation in Sec. III. Based on the traditional flux over population 
method we derive in Sec. IV the generalized expression for the Kramers' rate of escape 
from a metastable well. In Sec. V we numerically analyze the model and the bath modulated 
dynamics for the full potential and verify the theoretical rate with numerical simulation. The 
paper is concluded in Sec. VI. 

II. THE SYSTEM-RESERVOIR MODEL : THE RESERVOIR MODULATED BY 

AN EXTERNAL NOISE 

We consider a classical particle of mass M is linearly coupled to a heat bath of 
harmonic oscillators driven by an external noise. The total Hamiltonian is given by [|1^] 

H=^ + V{X) + ^ E + ^^^'(^^ - 9^^f] + H^nt . (1) 

In Eq.(|l|), X and p are the co-ordinate and momentum of the system particle; {qi,Pi) are the 
variables associated with the z-th oscillator and uJi and rrii are the corresponding frequency 
and mass, respectively. QiX measures the interaction between the particle and the bath. 
V{x) is the potential energy of the particle. Hint is assumed to be of the form 

1 ^ 

Hint = Ki Qi e{t) . (2) 

^ i=l 

The coupling function measures the strength of interaction and e{t) is the external noise 
which we assume to be stationary and Gaussian with zero mean, i.e., (e(t))e = and is 
characterized by an arbitrary correlation function as follows; 

(e(t)e(t'))e = 2D^(t-t') . 

Here (. . .)e implies the averaging over the external noise. 



We then eliminate the bath degrees of freedom in the usual way p6| , p!8|JT9| , |30[] to obtain 
the following generalized Langevin equation 
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X = f , 

^ = -^'1 dt' lit - t') v{t') + f{t) + 7r(t) (3) 
[ while constructing Eq.(^) we have set M and equal to unity ] where 

N 

lit) = m^l cos^it . (4) 

i=l 

f{t) is the internal forcing function generated through the coupling between the system and 
the heat bath and is given by 

N 

f{t) = ^gi {[qi{0) - gix{0)] coscj^t + t;i(0) cjj sincjjt} . (5) 

In Eq.(|^), 7r(t) is a fluctuating force term due to the external noise e{t) and is given by 

vr(t) = - /* ^{t - t') tit') dt' , (6) 
Jo 

where 

N 

(p{t) = Y^ Qi Ki uji sin Uit . (7) 

i=l 

The form of Eq.(^ therefore suggests that the system is driven by two forcing functions 
/(t) and vr(t). /(t) depends on the initial conditions of the bath oscillators for a fixed choice 
of the initial condition of the system degrees of freedom. To define the statistical properties 
of /(t), we assume that the initial distribution is one in which the bath is equilibrated at 
t = in the presence of the system but in the absence of the external noise agency such that 
{fit)) = and {fit) fit')) = ksT^it - t'). 

Now, at t = 0+, the external noise agency is switched on and the bath is modulated by 
e(t). The system is governed by Eq.@, where, apart from the internal noise fit), another 
fluctuating force nit) appears, that depends on the external noise e(t). Therefore, one can 
define an effective noise ^(t)(= fit) + irit)) whose correlation is given by 

((^(t)^(t'))) = ksT-fit - t') + 2D f dt' t dt" vit - t')vit' - t")^it' - t") , (8) 



5 



along with ((^(t))) = 0, where ((. . .)) means we have taken two averages independently. It 
should be emphasized that the above relation (|^) is not a fluctuation-dissipation relation 
due to the appearance of the external noise intensity. Rather it serves as a thermodynamic 
consistency condition. 

Let us now digress a little bit about vr(t). The statistical properties of 7r(t) are determined 
by the normal mode density of the bath frequencies, the coupling of the system with the 
bath, the coupling of the bath with the external noise and the external noise itself. Eq.@ 
is reminiscent of the familiar linear relation between the polarization and the external fleld 
where vr and e play the role of the former and the later, respectively. ip{t) can then be 
interpreted as a response function of the reservoir due to external noise e(t). The very 
structure of ^{t) suggests that this forcing function although originating from an external 
force is different from a direct driving force acting on the system. The distinction lies at the 
very nature of the bath characteristics (rather than system characteristics) as reflected in 
the relation and (0). 

With the coupling coefficients g{uj) = go/^/TcU and k,{uj) = .^TcUJKq, in the continuum 
limit PD| 7(t) and (f{t) reduce to the following forms 



= ^ I dcuV (to) cos ut (9) 



Tr 



and 

V'(^) = do 1^0 J duj 'D(uj) u! sincut . (10) 

where go and kq are constants and is the cutoff frequency of the oscillator. T>{uj) is the 
density of modes of the heat bath. 
From the above two relations, we obtain 

? = -Vit) ■ (11) 

Eq.(pTD is an important content of the present model. This expresses how the dissipative 
kernel 7(t) depends on the response function (f{t) of the medium due to external noise e{t) 



[ see Eq.(|^) ]. Such a relation for the open system can be anticipated in view of the fact 
that both the dissipation and the response function crucially depend on the properties of the 
reservoir especially on its density of modes and its coupling to the system and the external 
noise source. In what follows we shall be concerned with the consequences of this relation 
in terms of the Langevin description in the next section ( Eq.(|T^) ) and numerical analysis 
of the full model potential in Sec. V. 

III. THE FOKKER-PLANCK DESCRIPTION OF THE LINEARIZED MOTION : 
ASYMPTOTIC ANALYSIS OF THE FOKKER-PLANCK COEFFICIENTS 

We now consider the system to be a harmonically bound particle of unit mass and of 
frequency uq. Then because of Eq.(D the Langevin equation (H) becomes 

X = V , 

V = -u^x - t dt' -iit - t') v{t') + fit) - f dt' ^{t - t') t{t') (12) 
Jo Jo 

The Laplace transform of Eq. (|T^ allows us to write a formal solution for the displacement 
of the form 

x{t) = {{x{t))) + f dt' h{t - t') fit') - — r,u;2 m _ 

90 



go Jo 



-— r, f' dt' h^it - t') tit') , (13) 
9o Jo 



where we have made use of the relation ( [TT| ) explicitly. 
Here 

{{xit))) = x.it)m + hit)viO) (14) 

with x(0) and t'(O) being the initial position and initial velocity of the oscillator, respectively, 
which are nonrandom and 

xAt) = h-^l fhiT)dT] . (15) 
L Jo J 
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The kernel h{t) is the Laplace inversion of 

h{s) = -\ 2 (16) 

where, 7(5) = e~^* 7(t) dt, is the Laplace transform of the friction kernel 7(t), and 

hit) = . (17) 

The time derivative of Eq.(^) yields 

m = {{vm+ fdt' h^it-t') f{t')-^r,ul fdt' h{t-t') e{t') 
Jo go Jo 

-^rJ'dt'h3{t-t')e{t') (18) 
go Jo 



where 



{Mi)» = --'oM*) + i'(o)fti(«) , (19) 

= m = 'iM . (20) 

Next we calculate the variances. From the formal solution of x{t) and v{t), the explicit 
expressions for the variances are obtained which are given below; 

= 2 1'dh Kh) f ' dt, Hh) {f{ti)f{t2)) 

dh Kh) r dt, hit,) (e(ti)e(t2))e 

















{ 

Jo 











(21) 



<yi{t) = mt) - mt)))f)) 

= 2 fdu hit,) r dt, hit,) {fiti)fit,)) 

Jo Jo 
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(22) 



and 



criM = mt) - {{xmMt) - {{vit)))])) 



(23) 



where we have assumed that the noises f{t) and e{t) are symmetric with respect to the time 
argument and have made use the fact that f{t) and e{t) are uncorrelated. 

Due to the Gaussian property of the noises f{t) and e{t) and the hnearity of the Langevin 
equation (|12D, we see that the joint probabihty density p{x,v,t) of the oscillator must be 
Gaussian. The joint characteristic function associated with the density is 

1 



p(/i, p, t) = exp \ i{{x{t)))fi + i{{v{t)))p - ^ [al^{t)p'^ + 2al^{t)pp + al^{t)p^ 



(24) 



Using the method of characteristic function pT] , |52| and the above expression (^) we find the 
general Fokker-Planck equation associated with the probability density function p(x, f , t) for 
the process (|I^); 



where 



7(^) = -^lnT(t) , uJ^{t) = and 



(25) 



T(t) 



1 — UJ, 



lfdrh{T) 
Jo 



+ h\t) . 



The functions 0(t) and il){t) are defined by 

(pit) = ul{t)al^ + ^al^ + ]^al^ and ^it) = + ^{t)al^ + ulal^ - al^ 



(26) 



where the covariances are to be calculated for a particular given noise process. 

For the internal noise processes it had been shown earlier that for several models the 
various time dependent parameters ujQ{t), 7(t), etc. do exist asymptotically as t oo. The 
above consideration shows that h{t), hi{t), etc. do not depend on the nature of the noise 
but depend only on the relaxation ■y(t). 

We now discuss the asymptotic properties of 0(t) and ip{t), which in turn are dependent 
on the variances crl^it) and as t — » oo since they play a significant role in our further 

analysis that follows. 
From Eqs. (0) and (^2]) , we may write 

^lAt) = o^(t) + ol^\t) and aUt) = al^\t) + al^:\t) . 

where 'i' denotes the part corresponding to internal noise f{t) and 'e' corresponds to the 
external noise e{t). Since the average velocity of the oscillator is zero as t ^ oo we see from 
Eq.(p!9D that h{t) and hi{t) must be zero as t ^ oo. Also from Eq.(|14|) we observe that 
the function Xx{i) must decay to zero for long times. Hence, from Eq.(|T5|) we see that the 
stationary value of the integral of h{t) is I/liJq, i.e., 

/■oo 1 

/ h{t) dt = — . (27) 

Jo LUq 

Now, and (T2«(t) of Eqs.(|2l|) and (H) can be written in the form 

a^^it) = 2 f'dt, hit,) f " dt, h{t,) {f{h)f{h)) 

(28) 



C dr h{T) - h\t) - ^* dr h{T) 



and 



alf{t) = 2 fdt, hih) r dh hi{h) (/(ti)/(t2)) 
JO Jo 



1 - hiit) - uj'Qh\t) . (29) 



From the above two expressions [ Eqs.(p8|) and (|29| ) ] we see that 
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a2W(oo) = ^ and all}\oo) = ksT . (30) 

It is important to note that these stationary values are not related to the intensity and 
correlation time of the internal noise. 

We next consider the parts, and allj^\t), due to the presence of the external noise. 

The Laplace transform of Eq. ([TsD yields the expression 

x{s) - {{x{s))) = h{s)f{s) - ^r,u;fh{s)e{s) - ^T,s'h{sre{s) (31) 

where 

{{x{s))) = <-- otI 3^(0) + r^T 5 v{0) 

= |^-o;o'^|a;(0) + /i(s)t;(0) . (32) 

From the above equation (^) we can calculate the variance a^^ in the Laplace-transformed 
space which can be identified as the Laplace transform of Eq . (|2T|) . Thus, for the part o"^i^''(t) 
we observe that, cr^i'^^(s) contains terms like (^^TcUjQh{s)^ (e^(s))e. Since, we have assumed 
the stationarity of the noise e(t), we conclude that if C{0) exists [ where C(t—t') = {e(t)e(t'))e 
], then the stationary value of cr^^xK^) exists and becomes a constant that depends on the 
correlation time and the strength of the noise. Similar argument is also valid for a'^l^\t). 
Summarizing the above discussions we note that, 

(i) the internal noise-driven parts of crl^^t) and cr^„(t), i.e., a^^'^ and o"^^*^ approach the fixed 
values which are independent of the noise correlation and the intensity as t — oo, 

(ii) the external noise driven parts of variances also approach the constant values at the 
stationary {t oo) limit which are dependent on the strength and the correlation time of 
the noise. 

Hence we conclude, following the Ref.(33) and our preceding discussions that even in 
presence of an external noise the above terms do exist asymptotically and we write the 
steady state Fokker-Planck equation for the asymptotic values of the parameters as, 
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where, tUg, 7, 4>{oo), ip{oo), etc. are to be calculated from the general definition ( ]26|) for the 
steady state. As an explicit example we consider the case of a 5-correlated external noise 
and Ornstein-Uhlenbeck internal noise for which we provide the expressions for variances 
(j^^(t), and crlvit) and the relaxation function h(t) given in the Appendix-A. 

The general steady state solution of the above equation (Rl) is 



Pst{x,v) = -|exp 



2 — 2 

V LOqX 



+ 



2Do 2(Do + ^(oo))^ 



(34) 



where 



D, = ^ (35) 
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and Z is the normalization constant. The solution ([3^ ) can be verified by direct substitution. 
The distribution (|3^) is not an equilibrium distribution. This stationary distribution for the 
nonequilibrium open system plays the role of an equilibrium distribution of the closed system 
which may, however, be recovered in the absence of external noise term. 

IV. KRAMERS' ESCAPE RATE 

We now turn to the problem of decay of a metastable state. In Kramers approach the 
particle coordinate x corresponds to the reaction coordinate and its values at the minima of 
the potential V{x) denotes the reactant and product states. 

Linearizing the motion around barrier top at x = Xb the Langevin equation (0) can be 
written down as 



y = v , 



v = ujly- dt' 7(t - t') vit') + fit) + 7r(t) , (36) 







where, y = x — x^ and the barrier frequency u"^ is defined by 

1 



V{y) = H - ^ujly' ; ul>0 . (37) 



Correspondingly the motion of the particle is governed by the Fokker-Planck equation 
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dp 
di 



-vjf - u^l{t)yjf + %{t) — {vp) + + Mt)jri- , 

oy ov ov ov^ ovoy 



(38) 



where, the suffix '6' indicates that all the coefficients are to be calculated using the general 
definition (^) for the barrier top region. 

It is apparent from Eqs. (|33D and (|38D that since the dynamics is non-Markovian and 
the system is thermodynamically open one has to deal with the renormalized frequencies 
ujq and ujij near the bottom or top of the well, respectively. We make the ansatz that the 
nonequilibrium, steady state probability p^-, generating a nonvanishing diffusion current j, 
across the barrier is given by 

•2 V{X) 



Pb[X,V) 



exp 



2D, ^ D, + M^) 



(39) 



where 



76 



(40) 



V{x) is the renormalized linear potential as 

1 

2 
1 



V{x) = V{xo) + ^ujI{x - xof 



near the bottom 



V{x) = V{xb) — -f^lix — Xbf' , near the top 



(41) 



with ci}g, ul > 0. The unknown function $,{x,v) obeys the natural boundary condition that 
for X — > oo, ^{x,v) vanishes. 

The ansatz of the form ( pQ]) denoting the steady state distribution is motivated by the 
local analysis near the bottom and the top of the barrier in the Kramers' sense For 
a stationary nonequilibrium system, on the other hand, the relative population of the two 
regions, in general, depends on the global properties of the potential leading to an additional 
factor in the rate expression. Although because of the Kramers' type ansatz [jl| which is valid 
for the local analysis, such a consideration is outside the scope of the present treatment, we 
point out a distinctive feature in the ansatz (|39|) compared to Kramers' ansatz. While in 
the latter case one considers a complete factorization of the equilibrium part (Boltzmann) 
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and the dynamical part, the ansatz ( p9D incorporates the additional dynamical contribution 
through dissipation and strength of the noise into the exponential part. This modification of 
Kramers' ansatz (by dynamics) is due to nonequilibrium nature of the system. Thus unlike 
Kramers', the exponential factors in ( pQ]) and in the stationary distribution (|3^ ) which serves 
as a boundary condition are markedly different. Before carrying out global analysis in the 
present section our aim here is to understand the modification of the rate due to modulation 
of the bath driven by an external noise, within the perview of Kramers' type ansatz. The 
internal consistency of the treatment, however, can be checked by recovering the Kramers' 
result when the external noise is switched off. 

From equation (pHD, using ( pQ]) we obtain the equation for ^{y,v) in the steady state in 
the neighborhood of Xb, the equation 



Db / dy 



^by + ibv 



Db + ipb{oo) 

After making use of the appropriate transformations and boundary conditions for reduced 
distribution functions 0] we obtain the barrier crossing rate k given by 

-^0 



ujo_ Db ( A y/^ 



Db + V^b(oo) 



(43) 



where 



A 



A 



bb{oo) + ai)b{oo) 



with 



lb - V76 + 4t^6 r A = -7b - a 1 + 



2{Db + Moo)) i V ' - \ ' Db J - 

Here Eq is the activation energy, Eq = V{xb) — V{xo). Since the temperature due to internal 
thermal noise, the strength of the external noise and the damping constant are buried in the 
parameters Dq, Db, ipo, ipo and A the generalized expression look somewhat cumbersome. 
We point out that the subscripts '0' and '6' in D and ip refer to the well or barrier top 
region, respectively. Eg . (|43|) is one of the key results of this paper. We note here that 
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{Db + tlJh{oo))/kB in the exponential factor defines a new effective temperature characteristic 
of the steady state of the nonequihbrium open system and an effective transmission factor 
is contained in the prefactor controlhng the barrier crossing dynamics. As expected both 
are the functions of the external noise strength D and coupling of noise to the bath modes. 
The dependence of the rate on the parameters can be exposed explicitly once we consider 
the typical cases. 

V. NUMERICAL ANALYSIS OF THE MODEL AND THE BARRIER 

CROSSING RATE 

A. The model : bath modulation vs. direct driving 

We first explore the distinctive aspects of the external noise-driven-reservoir-modulated 
dynamics of the system ( in contact with the external noise driven bath ) in contrast to direct 
driving of the system by the external noise. This will help us to elucidate the special role of 
the reservoir response function in controlling the rate. To this end we first return to our basic 
equation of the model, Eq.(0) where (pit — t') as given by Eq.(|^ characterizes the response 
function. We solve this equation numerically using second order stochastic algorithm of Fox 
for the full potential V{x) = — and plot the results of computation of the inverse 
of the mean first passage time as a function of the damping constant Qq ( see Eq.(^) ) in 
Fig.(l). The density of modes of the bath has been assumed to be of the standard Lorentzian 
form, = ( ^2_|1^-^ ) with vanishing internal correlation time Tc. Typically for the 

other scaled parameters we have used fc^T = 0.1, the strength of the 5-correlated external 
noise, D = 0.1. The mean first passage times have been calculated by averaging over 10000 
trajectories. We first show the continuous curve in Fig.(l) which represents the case of 
direct driving by the external (5-correlated noise. What is immediately apparent is the loss 
of turnover of the rate as one allows the variation of the dissipation constant Qq from low 
to high friction regime. This is in sharp contrast to what is shown by the dotted curve in 
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Fig.(l) which depicts the situation of bath modulation with the response function if{t — t') 
determined by Eq.(|l^). Both the curves ( continuous and dotted ) representing the open 
system are compared to that for the standard Kramers' turnover, i.e., when the dynamics 
is calculated in absence of the external noise ( dashed curve ). As expected the rate in the 
driven system ( be it directly or through bath modulation ) is always higher that in the 
undriven one. Fig.(l) also shows that although in the low damping region the direct driving 
causes a much higher rate, it is, however, the bath modulation which becomes more effective 
in inducing activated barrier crossing in the high friction regime. It is thus interesting to 
note that the barrier crossing dynamics of the system in contact with an external noise 
modulated bath captures the essential turnover features of the Kramers' dynamics of the 
closed system. We therefore realize that although open, the nature of the response function 
of the reservoir as well as the thermodynamic consistency condition (H), make the open 
system feel like a closed system. 

In Fig. (2) we show the variation of inverse of the mean first passage time as a function 
of the strength of the external noise keeping all other parameter same as before but for 
(?Q = 2.0. It is apparent that the barrier crossing is more facilitated by modulating the bath 
than driving the system directly for higher values of external noise strength. 

B. The rate : Theoretical results vs. numerical simulation 

So far we have considered the full potential V{x) = — and vanishingly small 
correlation times for external and internal noises for numerical computation. We now turn 
to our basic theoretical result Eq.p3|) which is a generalization of Kramers' rate for bath 
modulated dynamics for intermediate to strong damping regime. To examine its validity 
we calculate the rate (|i3| ) as a function of the damping constant gl ( go is related to both 
7(t) and ip{t) in Eqs.(|^) and (p^Of ) respectively ) for several values of external noise strength 
D. The scaled barrier height Eq and ksT have been set to 0.25 and 0.1, respectively. The 
results are plotted in Fig. (3) ( continuous curve ) and compared to the rate, inverse of the 
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mean first passage time ( dotted curve ) calculated numerically using full potential for the 
set of parameter values ( as given in the earlier subsection ) with Eq.(^. In Figs. (4) and (5) 
we further compare our theoretical results ( continuous curve ) with numerical simulation 
( dotted curve ) for the variation of rate as a function of external noise strength and its 
correlation time, respectively for several values of damping constant. It is apparent from 
Figs. (3-5) that the theoretical and numerical results are in good agreement. 

VI. CONCLUSIONS 

Based on a system-reservoir microscopic model where the reservoir is modulated by an 
external, stationary and Gaussian noise with arbitrary decaying correlation function, we 
have numerically analyzed the model and generalized the Kramers' theory to calculate the 
steady state rate of escape from a metastable well. The main conclusions of this study are 
as follows; 

(i) We have shown that since the reservoir is driven by the external noise and the dis- 
sipative properties of the system depend on the reservoir, a simple connection between the 
dissipation and the response function of the medium due to the external noise can be estab- 
lished. This connection is important for realizing an effective temperature of the reservoir 
characterizing the stationary state of the thermodynamically open system, as well as an 
effective transmission factor controlling the rate. Both of these quantities depend on the 
strength and correlation time of the external noise. 

(ii) Many of the earlier treatments of the rate concern direct (phenomenological) driving 
of the system and did not emphasize the question of thermodynamic consistency in the con- 
text of open systems. The present theory being microscopic the fiuctuation-dissipation-like 
relation (P) remains an inbuilt characteristic of the model itself as an essential thermody- 
namic consistency condition. 

(iii) Based on numerical simulation of the full model potential we show that one can 
recover the turnover features of the Kramers' dynamics when the external noise modulates 
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the reservoir rather than the system directly. This recovery, we beheve, is an offshoot of the 
derived thermodynamic consistency condition (j^). 

(iv) Provided the long time limit of the moments for the stochastic processes pertaining 
to the external and internal noises characterized by arbitrary decaying correlation functions 
exist, the expression for generalized Kramers' rate of barrier crossing for the open system we 
derive here is fairly general. We have shown that it agrees reasonably well with numerical 
simulation using full potential for this model. 

The creation of a typical nonequilibrium open situation by modulating a bath with the 
help of an external noise is not an uncommon phenomenon in applications and industrial 
processing. The external agency generating noise does work on the bath by stirring, pump- 
ing, agitating, etc., to which the system dissipates internally. In the present treatment we 
are concerned with a nonequilibrium steady state which signifies a constant throughput of 
energy in contrast to thermal equilibrium defined by an constant temperature. We believe 
that these considerations are likely to be important in other related issues in nonequilibrium 
open systems and may serve as a basis for studying processes occurring within irreversibly 



driven environments P,P3[ and for thermal ratchet problems 
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APPENDIX A: CALCULATION OF VARIANCES 

We consider a particular case as an example where the external noise e{t) is 5-correlated 
and the internal noise is an Ornstein-Uhlenbeck process, i.e., 

{eit)e{t'))e = 2D6{t-t') 
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and 

(mm) = 

Consequently, from the fluctuation-dissipation relation we derive the dissipative kernel as, 

It should be noted that for Tc 0, the above noise process become (5-cor related. 
The Laplace transform of ^{t) as given above can be written as 

and subsequently, we have for Tc 7^ 0, 

s + a 



h{s) 



+ as"^ + bs + Co 
where 

1 , 2, ^0 1 ^0 
a — — , — uJq -\ and cq = — . 

''"c ''"c Tc 

We flnd that the inverse Laplace transform oih{s) reads 

h{t) = cie-^i* + cse""^^* sin(/3i + a) (Al) 
where the coefficients ci, C2, Ai, A2, /? and a are given by 

(A2a) 
(A2b) 

(A2c) 
(A2d) 

(A2e) 

1/3 

(A2f) 









A2 = 


l(A^B)^l, 












1 




2A2-Ai-d ' 






a(2A2 - Al) - A| - 




a - Al 




A = 


/ ab Cq 
-27 + y- 2 +^ 
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B 

C2 



a 
27 



Cl 



^ ab Cq 

2 1 ,5211/2 



1/3 



Q 



a = tan 



d-A2 



and 



2 1,2 



Q 



ah 



108 27 ^ 27 



6^4 



(A2g) 
(A2h) 

(A2i) 
(A2j) 



Here we note that for a physically allowed solution Ai, A2 must be positive. Since by Eq.(|T6D 
h{t) depends on the memory kernel 7(t) which is of decaying type and all the moments, in 
general, reach asymptotic constancy as shown in Sec. Ill, these quantities are positive (which 
depends on the correlation time Tc, the strength of the noise and other potential parameters) 
which may be checked (after some algebra) by considering the limiting cases such as Tc ^ 
and Tc large. 

Substituting Eq.( [AlD into the expressions for variances [external noise is 5-correlated] , 
namely into (^TJ) and (|2^) we have after some lengthy algebra 



2(e)/ 



where 



alf{t) = ksT[c,R + 



with 



Ai 



2 - CU^ C2i? + 



Cl 



Cl 



Ai 

2 - 2ujIc2R 



Ai 



+ e 



-Alt 



A1C1 + 



Ai , 



2c2e-^2t 
■Ai + /?2 



l-ulcR + ^^e-"^^' 
^1 



[A2 sin(/3t + a) + l3cos{i3t + a)] 



-2ciC2e-(^^+'^^)* sin(/?t + a) 
■ sm2(pt + a) 



2, ,2^2„-2A2i 



+ 



(Ai + P 



2\2 



(Ai + 



'2^2 



A^) 



(Ai + p^Y 



cie-^^^'sin\pt + a) 



(A3) 



R 



Al + P'- 



■(A2 sin a + /3 cos a) 



(A4) 
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and 



^l^:\t) = 2D (^^r,j [clU + At + 2ulAl)lA{t) 

-r2cM^t + A?(A^ - (5^) + ulil^l + - (5^)}Ic{t) 

-2c^/5A2(A2 - ^ ujl)lD{t) + 4c2/32A2/^(t) _ 4C1C2/3A2(A2 + Ujl)lF{t)] 
Here the /'s are defined by 



(A5) 



Similarly 



lA{t)= fe 
Jo 


-^^1* dt , 


(A6a) 


Jo 


-^"^'^ sm\pt + a) dt , 


(A6b) 


Jo 


-(^i+^2)*sin(/5t + a) dt , 


(A6c) 


lD{t)^ f e 
Jo 


-2^2*sin2(/3t + Q;) dt , 


(A6d) 


lE{t)= fe 
Jo 


-"^^^^ dt and 


(A6e) 


Irit)^ fe 
Jo 


-(^i+^2)*cos(/3i + Q;) dt . 


(A6f) 



<M-^lf{t) + al^:\t) 



where 



alf{t) = kBT - [(A? + u;o')c?e-2^^* + /J^c^e"^^^* 

-/^Asc^e-^^^* sin 2(/3t + a) + (A^ + a;^ - (3^)cle-^^''' sin\(3t + a) 
^g-(Ai+A2)t|2c,c2(a;^ + A1A2) sin(/3t + a) - 2Ai/3ciC2 cos(/3t + a)}] 



(A7) 



and 



[00 n 



?,Ai-PyP']lB{t) 



+2ciC2AiA2(u;o' + A^K' - 3/3' + A2)/c(i) 
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-2ciC2/3Ai(A? + u;1){ujI + SAj - f3')lF{t)] (A8) 



where, /'s are defined in Eq.( |A6a[ - [A61| ). The explicit expression for crl^it) can be derived 



from Eq.(^). In the hmit t ^ oo we calculate the stationary values of the variances. The 
variances cr^^(oo), cr^^(oo) and o"^^(oo) yield 0(oo) and ip^oo) and other relevant quantities. 
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FIGURES 

FIG. 1. Plot of barrier crossing rate, k vs. damping constant, g^. The solid and the dotted 
line correspond to direct driving of the system and bath modulation, respectively (D = 0.1). The 
dashed line corresponds to the thermodynamically closed system, i.e., the system without any 
external driving {D = 0.0). ksT = 0.1 is common for all the three curves, (units are arbitrary) 

FIG. 2. Plot of barrier crossing rate, k vs. external noise strenth, D for a constant g^. The 
curve (a) represents the results for bath modulation while the curve (b) is the result for direct 
additive driving, (units are arbitrary) 

FIG. 3. Plot of barrier crossing rate, k vs. damping constant, for different external noise 
strengths D. The solid lines correspond to theoretical result ( Eq.(52) ) and the dotted curves are 
due to simulation, (a) D = 0.15, (b) D = 0.10 and (c) D = 0.05. (units are arbitrary) 

FIG. 4. Plot of barrier crossing rate, k vs. external noise strength, D for different values of 
5q. The solid and the dotted lines are same as in Fig. (3). (a) g^ = 2.0 and (b) gg = 3.0 . (units 
are arbitrary) 

FIG. 5. Plot of barrier crossing rate, k vs. correlation time of the external noise, Tg for different 
values of external noise strengths, D. The solid and the dotted lines are same as in Fig. (3). (a) D 
= 0.5, (b) D = 1.0 and (c) D = 1.5 . (units are arbitrary) 
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